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^ ' Abstract 

' The design of sensor networks capable of reaching a consensus on a globally optimal decision test, 

f>f~j ' without the need for a fusion center, is a problem that has received considerable attention in the 
last years. Many consensus algorithms have been proposed, with convergence conditions depending 

' on the graph describing the interaction among the nodes. In most works, the graph is undirected 

■ and there are no propagation delays. Only recently, the analysis has been extended to consensus 

^ ■ algorithms incorporating propagation delays. In this work, we propose a consensus algorithm able 
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to converge to a globally optimal decision statistic, using a wideband wireless network, governed 
by a fairly simple MAC mechanism, where each link is a multipath, frequency-selective, channel. 
The main contribution of the paper is to derive necessary and sufficient conditions on the network 



> 

, topology and sufficient conditions on the channel transfer functions guaranteeing the exponential 

convergence of the consensus algorithm to a globally optimal decision value, for any bounded delay 
condition. 
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; 1 Introduction 
> . 

^ ■ Distributed algorithms for achieving consensus in wireless sensor networks, without the need for a 
' fusion center, have been the subject of many recent works. Two excellent tutorials on the subject are 
[U [2] (see also references therein) . The conditions for achieving a consensus over a globally optimal 
decision test ultimately depend on the properties of the graph modeling the interaction among the 
nodes. Most works consider undirected graphs and neglect propagation delays. There are only a few 
works that study the impact of delays in consensus-achieving algorithms, namely [3] — [B], focusing on 
time-continuous systems, and [7] — [9], dealing with discrete-time systems. Among these works, it is 
useful to distinguish between consensus algorithms, H] — [3], where the states of all the sensors converge 
to a prescribed function (typically the average) of the sensors' initial values, and agreement algorithms, 
[Ij — [9], typically used for coordinating the motion of sets of vehicles, where the states of the nodes 
converge to a common value, but this value is not a specified function of the initial values. A recent 
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work proposed a randomized gossip algorithm [TO] to achieve distributed consensus, with a simple 
interaction mechanism, where each node interacts with one node at the time, in a randomized fashion. 

In this work, we are interested in distributed consensus algorithms where the consensus coincides 
with a globally optimal decision statistic. Our goal is to derive the conditions on the channels between 
each pair of nodes, guaranteeing that each sensor will eventually converge to the globally optimal 
decision statistic, in a totally distributed manner, i.e. without requiring the presence of a fusion 
center. In [H [3], the authors provided necessary and sufficient conditions for the convergence of a 
linear consensus protocol, in the case of a common time-invariant delay value for all the links, i.e., 
Tij = T Vi ^ j, and assuming symmetric channels among the nodes (modeled as a undirected graph). 
Under these assumptions, the average consensus in [ll [3] is reached if and only if the common delay 
r is smaller than a topology-dependent value. However, the assumptions of homogeneous delays 
and nonreciprocal channels are not appropriate for describing the propagation in a common network 
deploying scenario, where the delays depend on traveled distances and the communication channels 
may be asymmetric. In [TT], we generalized the consensus algorithms to networks with inhomogeneous 
delays and asymmetric flat-fading channels. In this correspondence, we extend our previous work 
to the more general case where each link is modeled as a multipath channel. We assume baseband 
communications, motivated by the use of impulse radio technologies. The main contributions of this 
paper are the following: i) We provide necessary and sufficient conditions on the network topology 
and sufficient conditions on the transfer function of each channel ensuring global convergence to the 
optimal decision test, for any set of finite propagation delays; ii) We prove that the convergence is 
exponential, with convergence rate depending, in general, on the channel parameters and propagation 
delays; iii) We show how to reach a distributed consensus coinciding with the globally optimal decision 
statistics, achievable by a centralized system having error-free access to all the nodes measurements 
and observation parameters, without the need of estimating neither the channel coefficients nor the 
delays. 

2 How to Achieve Consensus on a Globally Optimal Decision Test 
in a Decentralized Way 

Let us consider a set of sensors, each measuring a scalar parameter yi, i = 1,. . . ,N. The goal of 
the network is to compute a sufficient statistic of the measured data expressible as 



where {cj} are positive coefficients and {gi} and h are arbitrary (possibly nonlinear) real functions on 
M, i.e., gfj, /i : M 1-^ M. Even though the class of functions expressible as in ([T]) is not the most general 
one, it does include many cases of practical interest, like, e.g., best linear unbiased estimation or ML 
estimation under linear signal models, multiple hypothesis testing, detection of Gaussian processes 
in Gaussian noise, computation of maximum, minimum, geometric mean or the histograms of the 
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gathered data |1H I13j. In this paper, we consider only the scalar observation case, but the extension 
of ([1]) to the vector case is straightforward, along the same guidelines of pT} [T2]. 

To compute functions in the form ([1]) in a distributed way, we consider a linear interaction model 
among the nodes, and we generalize the approach of |1H [T2] to a network where the channel be- 
tween each pair of nodes is a multipath channel, with, in general, asymmetric channel coefficients and 
geometry-dependent delays. In each node there is a dynamical system whose state Xi{t;y) evolves 
according to the following linear differential equation 



where y = {yi}fLi is the set of measurements; gi{yi) is a function of the local measurement, whose 
form depends on the specific decision test; Cj is a positive coefficient that is chosen in order to achieve 
the desired consensus, as in ([1]); K is a positive coefficient controlling the convergence rate; a« and 
rf^ are the amplitude and the delay associated to the l-ih path of the channel between nodes i and 
j; Mi = {j = 1, . . . , N : 3 afj 7^ 0, / = 1, . . . , L} denotes the set of neighbors of node i, i.e., the nodes 
that send signals to node i. It is worth noticing that the state function of, let us say, node i depends, 
directly, only on the measurement yi taken by the node itself and only indirectly on the measurements 
gathered by the other nodes. In other words, even though the state Xi{t; y) gets to depend, eventually, 
on all the measurements, through the interaction with the other nodes, each node needs to know only 
its own measurement. 

The channel through which node r receives the signal from node g is a multipath channel with 
transfer function Hrq{juj) = X^^i argC"-^'^'^^'' , for all r ^ q. We assume that the channel coefhcients 
are sufficiently slowly varying to be considered constant for the time interval necessary for the network 
to converge, within a prescribed accuracy. In Section 3, we will show that the convergence of ([2|) is 
exponential and we will derive a bound for the convergence rate. Knowing this rate, our method is 
applicable for those channels whose coherence time is sufficiently greater than the convergence time. 
We are interested in baseband communications, motivated from a possible implementation of the 
radio interface allowing for the interaction described by ([2]) with an impulse radio using pulse-position 
modulation (IR-PPM), where the position of the pulse transmitted by node i is proportional to the 
state of node Xi{t). In general, we allow the channels to be asymmetric, i.e., a^q may be different from 
Ogr (and thus Hrq{juj) 7^ Hqj.{juj)). We also assume, realistically, that the maximum delay is bounded, 
with maximum value r = max^^g^; Trq ■ Because of the delays, the state evolution ([2]) for, let us say, 
t > 0, is uniquely defined provided that the initial state variables Xi{t;y) are specified in the interval 
from — T to 0, i.e., Xi{i9] y) = </>j(i?), for all i = 1, . . . , A^, and -d E [— r, 0]. 

Some important comments about the interaction mechanism ([2]) are appropriate. Distributed 
consensus algorithms have a clear advantage with respect to centralized systems, as they are less prone 
to congestion events or failures of some of the nodes. They are also inherently scalable. However, 
as opposed to centralized systems, they typically require an iterative mechanism to converge to the 




* j£N,. 1=1 
x^i^;y)=U^), i9g[-t, 0] 
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desired decision test. In most available works on distributed consensus, it is tacitly assumed that 
each node is able to receive the signals sent by its neighbors separately. This, of course, requires 
a proper medium access control (MAC) mechanism to avoid collisions. But, when combined with 
the iterative nature of distributed consensus algorithms, a collision avoidance MAC protocol may 
become rather complicated. Even the simple randomized gossip algorithm of jlDj requires some form 
of MAC to avoid collisions. Unfortunately, enforcing a MAC control goes against the requirement of 
simplicity and scalability, which are some of the major motivations underlying the use of distributed 
consensus algorithms. Conversely, we are interested in distributed consensus mechanisms where all 
nodes transmit over a common shared physical channel and there are no collision avoidance or resolution 
mechanisms whatsoever, so that each node receives a linear combination of the signals transmitted by 
the other nodes, possibly through a multipath propagation channel. This motivates the interaction 
model expressed by from which it turns out that each node does not need to resolve the received 
signals to be able to update its own state function. In this correspondence, we do not study the radio 
interface allowing for the node interaction given by ([2]). Nevertheless, some preliminary studies, see 
e.g., [I3l[20l[21] suggest that impulse radios with pulse position modulation or distributed phase- lock 
circuits are possible candidates for implementing where the state values are exchanged through 
pulse position modulation or phase modulation, respectively. 

However, the advantages of distributed consensus based algorithms as described above come at 
the price of a penalty: the final consensus is reached through an iterative procedure that consumes 
time and energy. The overall energy necessary to achieve the final decision is the sum of the powers 
transmitted by each sensor multiplied by the convergence time (in the next section, we will give an 
upper bound of such a value). On one hand, to save energy, we would like to use the minimum 
transmit power that ensures network connectivity. But a small transmit power has an effect on the 
network topology, as it leads to a reduced number of links and, as a consequence, to a small algebraic 
connectivity. Hence, a small individual transmit power implies a long convergence time. Conversely, 
to reduce the convergence time, the network should have a high connectivity, but this requires a large 
transmit power. It is then intuitive to expect an optimal trade-off. This trade-off has been studied 
in [22], where we remand to the interested reader. The focus of this correspondence is on finding the 
conditions on the channel parameters that guarantee the convergence of ([2]) to the desired consensus 
value. 

Consensus on the state derivative. Differently from most papers dealing with average consensus 
problems [I] — [3] , [1] — [9] , we adopt here the alternative definition of consensus already introduced in 
our previous works — [T3]: We define the consensus (or network synchronization) with respect to 
the state derivative, rather than to the state. 

Definition 1 Given the dynamical system in we say that a solution {x*{t]y)} of ^ is a synchro- 
nized state of the system, if x*{t; y) = a*{y), Vi = 1,2, . . . , N . The system (0) is said to globally syn- 
chronize if there exists a synchronized state a*{y), and all the state derivatives asymptotically converge 
to this common value, for any given set of initial conditions {0j}, i.e., lim^H^oo \xi{t;y) — a*{y)\ = 0, 
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Vi = 1,2, ...,A'", where {xi{t;y)} is a solution to The synchronized state is said to be globally 
asymptotically stable if the system globally synchronizes, in the sense specified above. 

Observe that, according to Definition [H if there exists a globally asymptotically stable synchronized 
state, then it must necessarily be unique (in the derivative). One of the reasons to introduce this 
definition of consensus, as opposed to the consensus on the state [l] — [3], [4] — [9], is that, as will be 
shown in the next section, the convergence on the state derivative is not affected by the presence 
of propagation delays. One more reason is that, in the presence of coupling noise, state-convergent 
algorithms give rise to a noise with diverging variance , whereas the algorithm converging on the state 
derivative exhibits a finite variance [12\ [T3] . 



3 Necessary and Sufficient Conditions for Achieving Consensus 

To derive our main results, we rely on some basic notions of directed graph (digraph) theory, as briefly 
recalled next. More details are given in }1H Appendix A]. A digraph ^ is defined as ={'^, where 
is the set of vertices and C x is the set of edges, with the convention that = {vi,Vj) G co 
if there exists an edge from Vj to Vi, i.e., the information flows from Vj to Vi. A digraph is weighted 
if a positive weight, denoted by Ojj, is associated with each edge Cij. The in-degree of a vertex is 
defined as the sum of the weights of all its incoming edges. The out-degree is similarly defined. The 
Laplacian matrix L = L(?^) of the digraph associated to system ([2]) is L = D — A, where D is the 
diagonal matrix of vertex in-degrees and A is the adjacency matrix. For reasons that will be clarified 
within the proof of next theorem, the above matrices are built as follows: [Dju = YljeMt X^ti '^fj ^"^^ 
[A\ij = X^^x '^^ij- ^ digraph is a directed tree if it has N vertices and — 1 edges and there exists a 
root vertex (i.e., a zero in-degree vertex) with directed paths to all other vertices. A directed tree is a 
spanning directed tree of a digraph ^ if it has the same vertices of A digraph is Strongly Connected 
(SC) if, for every pair of nodes Vi and Vj , there exists a directed path from vi to Vj and viceversa. A 
digraph is Quasi- Strongly Connected (QSC) if, for every pair of nodes vi and Vj , there exists a node 
r that can reach both Vi and Vj by a directed path. The fundamental result of this paper is stated in 
the following. 

Theorem 1 Let L be the Laplacian matrix associated to the digraph ^ ={^, of system (0), and 
let'y = [7;^, . . . , 7^]^ be the left eigenvector of L corresponding to the zero eigenvalue, i.e., 7"^L = 0^. 
Given system assume that the following conditions are satisfied: 
al) The coupling gain K and the coefficients {q} are positive; 

a2) The propagation delays {^fj^ ^'"e finite, the coefficients {O'^ij} O'f^ '^ea/ and the channel transfer 
functions {Hrq{juj)} are such that 

Hr,iO)>0, yq,r^q, and ^^^^l}^^^^ < 1 , Va;GM,Vr/g; (3) 

a3) The initial conditions are taken in the set of continuously differentiable and bounded functions 
mapping the interval [— r, 0] to M^. 
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Then, system (0) globally synchronizes, for any set of propagation delays, if and only if the digraph 
^ is QSC. The synchronized state is 



where 'j^ > if and only if node i can reach all the other nodes of the digraph by a directed path, 
otherwise 7^ = 0. The convergence is exponential, with asymptotic convergence rate arbitrarily close 
to r = — minj{|Re{si}| : p{si) = and Si 7^ 0}, where p{s) is the characteristic function associated to 
system (0j (see (W(l in the Appendix). 



Remark 1 - Robustness against multipath channels: Theorem [T] shows that, differently from 
classical linear consensus protocols [HE], the proposed algorithm is robust against propagation delays, 
since its convergence condition is not affected by the delays. Moreover, the proposed approach is valid 
for frequency- selective and asymmetric channels. The only significant constraint is that the channel 
coefficients are real and that, according to a2), their summation, over each channel, has to be a positive 
quantity. This implies a sort of implicit coherent combination, conceptually similar to the type-based 
approach of |19j . even though the work of [19j was aimed at studying the multiple access for sensor 
networks with a fusion center, whereas our scheme does not need a fusion center. The additional 
constraints on the channel transfer functions, i.e., '^ZqeNr W'rq{j^)\ < YlqeMr ^rq{0), for all w G M, is 
certainly valid if the channels are low-pass filters with maximum gain in w = 0. But this is only a 
sufficient condition and we will later report some numerical results showing that if this condition is 
not satisfied, the method can still converge. Finally, given the convergence rate r, we know for which 
class of channels the method is applicable: the channels whose coherence time is sufficiently greater 
than 1/r. 

Remark 2 - Effect of network topology: According to Theorem [H a global consensus is possible if 
and only if there exists at least one node (the root node of the spanning directed tree of the digraph) 
that can reach all the other nodes by a directed path. If no such a node exists, the information gathered 
by each sensor has no way to propagate through the whole network and thus a global consensus cannot 
be reached. Moreover, the only nodes contributing to the final consensus value are the ones having a 
directed path linking them to all the other nodes [see ^] . As a consequence, the final decision depends 
on the measurements gathered by all the nodes if and only if the network is strongly connected. When 
the digraph is not QSC, system ([2]) may still converge, but it forms separated clusters of consensus, as 
proved in [TTl [13] . 

Remark 3 - Unbiased decisions without estimating the channel parameters: The closed 
form expression of the synchronized state given in shows a dependence of the final consensus on 
the network topology and propagation parameters. This implies that the final consensus value dH), 
in general, does not coincide with the desired decision statistics as given in ([T]), except for the trivial 
case of fiat-fading channels with zero delays and balanced (and thus strongly connected) digraph. 





Proof. See the Appendix. 
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Nevertheless, in the following we provide a method to get an unbiased estimate, without having to get 
any preliminary estimation of the channel parameters, i.e. incorporating, only in the case of 

unbalanced networks, a decentralized estimation of the topology dependent coefficients 7^. 

The bias due to the propagation delays and path amplitudes can be removed using the following 
two-step algorithm. We let system ([2]) to evolve twice: The first time, the system evolves according to 
([2]) and we denote by a*{y) the synchronized state, as in the second time, we set gi{yi) = 1 in ([2]), 
for all i, and the system is let to evolve again, denoting the new synchronized state by a*(l). Taking 
the ratio a*{y)/a*{l) = {Yli^=i li'^idiiyi)) / Cl2iLi 7iCi)> we obtain the same consensus value that would 
have been achieved in the absence of multipath propagation. 

If the network is strongly connected and balanced, 7j = 1, Vi and then the compensated consen- 
sus coincides with the desired value ([T]). If the network is unbalanced, the compensated consensus 
a*{y)/a*{l) does not depend on the multipath coefficients, but it is still biased, with a bias dependent 
on 7, i.e., on the network topology. This residual dependence can be eliminated in a decentralized 
way if each node is able to estimate its own 7^. In fact, in such a case, a*{y)/a*{l) can be made to 
coincide with the desired expression in ([T|) by simply replacing each Cj in ([2]) with Cj/7j, for all i such 
that 7j / (suppose that there are Nr of such nodes, w.l.o.g.). Interestingly, the estimate of each 7, 
can also be obtained in a decentralized way, using the following procedure. At the beginning, every 
node sets gi{yi) = 1 and the network is let to evolve. The final consensus value will be, in this case 
a*(l). Then, the network is let to evolve Nr times, according to the following protocol. At step i, 
with i = 1, . . . Nr, node i sets gi{yi) = 1, while all the other nodes set gkivk) = for all k 7^ i; all 
nodes are then let to evolve according to ([2|); let us denote by a*{ei) the final consensus value, where 
Cj is the canonical vector having all zeros, except the i-th component, equal to one. Each node is now 
able to take the ratio a*(ej)/a*(l), which coincides with the ratio 7, := 7j/ 7^^.. Thus, after Nr + 1 
steps, every node knows its own (normalized) 7j and it may then use it in the subsequent run of the 
consensus algorithm, setting q = c/7j, to achieve a topology independent estimate. Observe that, 
since the eigenvector 7 does not depend on the observations {yj}, the proposed algorithm to estimate 
7 is required to be performed only once every channel coherence period. In summary, the effects of 
both delays and channel coefficients can be eliminated from the final consensus value, even if at the 
price of a slight increase of complexity and the need for some coordination among the nodes. 

4 Numerical Results and Conclusion 

As a numerical example, in the top row of Fig. [H we report two examples of topologies: the left graph 
is SC, whereas the right graph is QSC. For the QSC digraph in the figure, we sketch its decomposition 
in Strongly Connected Components (SCC), whose root is denoted by RSCC0 The behavior of the 
state derivatives versus the iteration index, in the two cases, is illustrated on the bottom row of the 

^ A SCC of a digraph is a maximal subgraph which is also SC, meaning that there is no larger SC subgraph containing 
the nodes of the considered component. A RSCC is a SCC containing all nodes that can reach all the other nodes in the 
digraph by a directed path [llj Appendix A] . 
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figurqj. The edges shown in both graphs show the active hnk. Each hnk is modeled as an FIR fiher, 
modehng the multipath fading. Each filter has maximum length L = 5 and the coefficients have been 
generated as afj = {A + 'w{i,j, I)) e~'''^^'^°, Z = 0, . . . , L — 1, where the constant A = 1 represents a 
deterministic component, whereas w{i,j,l) are i.i.d. random Gaussian variables with zero mean and 
standard deviation cr„ = 0.5, modeling the fading. Observe that, using this setting, some channel 
coefficients are also negative. The exponential models the attenuation as a function of distance and tq 
represents the delay spread; T is the sampling time. The delays rfj on each link have been modeled as 
rfj = dij/c+ {I — 1)T, where dij is the distance between nodes i and j and c is the speed of light. The 
dimension of the network has been computed in order to make the maximum delay T^ax = dmax/c 
much larger than the sampling time T. In particular, we chose the parameters so that Tmax = 30 T, in 
order to test the algorithm under a severe propagation delay. The constant lines with arrows reported 
in the bottom row of Fig. [1] represent the theoretical value, as given by (j3]). We can verify that the 
simulation curves tend to approach the theoretical values for both SC and QSC topologies, as predicted 
by the theory. It is worth mentioning that, in both cases, we used channels that respect the condition 
Hrq{0) > 0,Vr, g, but do not necessarily respect the condition '^q^j\/^\Hrq{juj)\ < l^gg^^ -ffrij(O). 
Nonetheless, the simulation results are still in good agreement with our theoretical findings. There is 
no contrast with the theory because the condition J2qeAfr \^rq{j^)\ < J2qeMr ^rqi^), Vr 7^ q, is only a 
sufficient condition. 

The estimates reported in Fig. [1] show a good agreement between theory and simulation, but 
the final result does not coincide with the theoretical optimal value, because of the bias induced by 
the multipath coefficients and delays. However, as suggested at the end of the previous section, it is 

(k) 

possible to remove the bias, without having to estimate neither the channel amplitudes J nor the 
delays r} • . As an example of this compensation technique, in Fig. [21 we report the running state 
derivative Xi{t;y) (solid lines) and the compensated estimate Xi{t;y)/xi{t;l) (dotted line), together 
with the theoretical limits (constant lines) achievable without compensation (triangle marks) and with 
compensation (circle marks). The last value coincides with the globally optimal estimate. The results 
shown in Fig. [2] have been achieved with multipath fading channels of length L = 11, under the 
same fading model used in the previous example. Fig. [2] shows that, as predicted by the theory, the 
consensus algorithm with compensation is able to reach the globally optimal estimate, without the 

need of estimating the channel coefficients. 

In summary, in this work we have derived the conditions allowing a distributed consensus mecha- 
nism to reach globally optimal decision statistics, in the presence of multipath propagation in the link 
between each pair of nodes. The method is valid for real (baseband) channels and it requires that the 
summation of the channel coefficients over each link is strictly positive. Thanks to the closed form 
expression derived in this paper, we have also shown how to get unbiased, globally optimal estimates, 
without the need to resolve the signals received from different nodes (thus allowing for a very simple 
MAC mechanism), nor to estimate the channel coefficients. A crucial investigation, motivated from 

^Clearly, the simulations have been performed on the discretized version of ([2]); in such a case, given the samphng time 
T, there is a maximum value of K guaranteeing the convergence of KT must be sufBciently smaller than 1/ deg™'''', 
where degj^"" is the maximum in-degree of the graph Laplacian (cf. [Ill Appendix A]). 
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Figure 1: Example of global consensus in a SC and QSC network. 

this work, is the design of the most appropriate radio interface allowing for the internode interaction 
enabhng the distributed consensus. 



5 Appendix 

Part of the proof of the theorem is based on the same approach we foUowed in |1H Appendix C]. Thus, 
in the following we make use of the general results of [llj Appendix C] and we focus only on the specific 
aspects of the model used in this paper. We start studying the existence of a synchronized state in 
the form Then, we prove that such a state is also globally asymptotically stable (cf. Definition 
1). Throughout the proof, we assume that conditions al)-a4) are satisfied and that the digraph ^ 
associated to ([2]) is QSC. In the following, for the sake of notation simplicity, we drop the dependence 
of the state function from the observation, as this dependence does not play any role in our proof. 

Existence of a synchronized state: The set of delayed differential equations ([2]) admits a solution 
in the form 

x*{t) = at + xlo, i = l,...,N, (5) 

where a € M and {a^^^oi ^ "-"^ coefficients that depend in general on the system parameters and 
on the initial conditions, if and only if satisfies ([2]), i.e., if and only if there exist a and {x*q} 
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Figure 2: Uncompensated running estimate (solid lines), i.e. Xi{t;y) and compensated running estimate (dotted lines), 
i.e. Xi{t; y)/xi{t\ 1). 



such that the following system of linear equations is feasible: 

K 



K 



j,0 -^1,0 



Vi = 1, . . . , A^, where 



Ai(a) = g^ivi) - a ( 1 + ^ ^ .^^^ 



^ (0 ,(0 

^3 «i 



1=1 



(6) 



(7) 



Introducing the weighted Laplacian L = L(l^) associated to 1^, system ([6]) can be equivalently rewritten 
in vector form as 

i^LxS = DeA(Q), (8) 



where = [xf g, • • • jX^g]"^, Dc = diag(ci, . . . ,CAr), and A 



a 



[Ai(a),...,A^(a)]^, with Ai(a) 



defined in ([7]). Observe that, because of a2) and the quasi-strong connectivity of Sf, the graph 
Laplacian L has the following properties: i) rank(L) = — 1; ii) A/'(L) = 1Z{\]^) ; and iii) Mijj^) = 
7^(7) , where M{-) and 7^(-) denote the (right) null-space and the range space operators, respectively, 
and 7 is a left eigenvector of L corresponding to the (simple) zero eigenvalue of L, i.e., 7"^L = 0"^. It 
follows from i)-iii) that, for any given q, ([8]) admits a solution if and only if DqA (a) € 7^(L). Using 
again properties i)-iii), we have: DcA(a) G '^(L) 7'^DcA(a) = 0. It is easy to check that the 
value of Q that satisfies the latter condition is a = a*, with a* defined in Hence, if q = a*, the 
synchronized state in the desired form ([5]) is a solution to ([2]), for any given set of {r^ij}, {Oi}, {cj}; 
{a^^} and K > Q. 

The structure of the left eigenvector 7 associated to the zero eigenvalue of L as given in the theorem 
follows from [TTl Lemma 4] . 

Setting a = a*, system 1^ admits 00^ solutions, given by Xg = -^L''DcA(q*) -|- TZ{1n) — xq -|- 
7^ (Iat) , where 

xo ^ L«DeA(a*)/i^, (9) 
Ai{a*) is obtained by ([7]) setting a = a* and is the generalized inverse of the Laplacian L. 

Global asymptotic stability: We prove now that the synchronized state of system ([2|) is globally 
asymptotically stable. To this end, we use the following intermediate results. 
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Let C+ = {s G C : Re{s} > 0}, C_ = {s € C : Re{s} < 0}, and C+ be the closure of C+, i.e^ 
C+ = {s G C : Re{s} > 0}. Denoting by H'^^"'- the set of n x m matrices whose entries are analytiqfl 
and bounded functions in C+, let us introduce the degree matrix A > (where ">"has to be intended 
component-wise) and the complex matrix H(s) £ C^^^, defined respectively as 

A A I 0' if ^ = i) 

A = diag(A;idegin(i^i),...,A;Ardegin(wAr)), [n{s)\- = I m (10) 

( h Ei=i % e , if « / j, 

where deg^^{vi) = J2jeJ\fi Yl[=i'^fj > is the in-degree of node vi and ki = K/ci > 0. Observe that 
H(s) G 

Lemma 1 Consider the following linear functional differential equation: 

Xi{'d) = M'd), i9e[-T, 0]. 

and assume that the following conditions are satisfied: 

hi. The initial value functions cj) are taken in the set of continuously difjerentiable functions 
that are bounded in the nor — sup_^<^<;Q 11^ , and the solutions x(t) with initial 

functions cf) are bounded; 

b2. The characteristic equation associated to flT]] 

= det(sl + A-H(s)) = 0, (12) 

with A and H(s) defined in fW\). has all roots {sr}r G C_, with at most one simple root at 
s = O0 

Then, system ill]) is marginally stable, i.e., \/cf) € C"^ and Re{si} < c < 0, there exist ti and a, 
with to < ti < +00 and < a < +00, independent of cf), and a vector x°°, with ||x°°j| < +00, such that 

||x(t) -x°°|| < a|0|,e^*, yt>ti. (13) 

Proof. Because of space limitation, we omit the proof that can be obtained following the same steps 
of the proof in [11^ Lemma 5] , after observing that system (jlip can be rewritten in the canonical form 
of dSl Ch. 6, Eq. (6.3.2)], [T71 Ch. 3, Eq. (3.1)]. ■ 

Lemma 2 ( [15L Theorem 2.21) Let H(s) G y^^xN p(H(s)) denote the spectral radius o/H(s). 
Then, p{ii{s)) is a subharmonic^ bounded (above) function on C+. □ 



complex function is said to be analytic (or holomorphic) on a region !> C C if it is complex differentiable at every 
point in V, i.e., for any zq £ "D the function satisfies the Cauchy-Riemann equations and has continuous first partial 
derivatives in the neighborhood of 20 (see, e.g., 1141 Theorem 11.2]). 

*We used, without loss of generality, as vector norm in the infinity norm IMI;,^ , defined as ||x||^ = max^ \xi\. Of 
course, the same conclusions can be obtained using any other norm. 

^We assume, w.l.o.g., that the roots {sr} are arranged in nonincreasing order with respect to the real part, i.e., 

= Re{so} > Re{si} > Re{s2} > .... 

^See, e.g., 1151 . for the definition of subharmonic function. 
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We are ready to prove the global asymptotic stability of the synchronized state of ([2]) . Applying the 
following change of variables: ^'i(t) — Xi{t) — [a*t + Xi^), for alH = 1, . . . , A^, where a* and {xj^o} are 
defined in (HD and dHJ, respectively, and using the original system ([2]) can be equivalently rewritten 
in terms of {^i{t)}i as 

= ^.^^^ 41 (^.(^ - 4) - ^«(*)) ' * ^ 0' (14) 

with = (/>j(i?) = </'i('!9) — a*^ — Xifi for ?? € [— r, 0], where {i;^)^} are the initial value functions of 

the original system ([2]). 

It follows from (jl4p that the synchronized state of system ([2]), as given in ^ (with a = a*), is 
globally asymptotically stable (according to Definition [T|) if system (|14p is marginally stable. According 
to Lemma [H the marginal stability of system (jl4p is guaranteed if: bl) the trajectories {^i{t)} are 
bounded for all t > 0, given cj) £ C^; b2) the characteristic equation (jl2p associated to (|14p . has all 
roots in C_, with at most one simple root at s = 0. 

Following the same steps as in |1H Appendix C], one can prove that, under al)-a3), all the solutions 
{^i{t)} to (jl4p . with initial conditions in C^, are uniformly bounded, as required by assumption bl) in 
Lemma[T] Because of space limitation we omit the details. We study instead the characteristic equation 
(jl2p . and prove that, under al)-a3) and the quasi-strong connectivity of the digraph, assumption b2) 
of Lemma [1] is satisfied. 

First of all, observe that, since A — H(0) = i^DcL, we have 

p(0) = det (A - H(0)) = det (iTDc) det (L) = 0, (15) 

where the last equality in (|15p is due to rank(L) = A — 1. It follows from (jlSp that p{s) has a root in 
s = 0, corresponding to the zero eigenvalue of the Laplacian L (recall that det (A'Dc) ^ 0). Since the 
digraph is assumed to be QSC, according to [11, Corollary 2], such a root is simple. Thus, to complete 
the proof, we need to show that p{s) does not have any solution in C-|-\{0}, i.e., 

det(sI+A-H(s)) /O, Vs G C+\{0}. (16) 

Since sI+A is nonsingular in C+\{0} [recall that, under al)-a2), A > 0, with at least one positive 
diagonal entry], (fTUj) is equivalent to 

det (l-(sI+A)"^H(s)) 7^ 0, Vs G C+\{0}, (17) 

which leads to the following sufficient condition for (fT6]) : 

pis)^p ((sI+A)-i H(s)) < 1, V.S G C+\{0}. (18) 

Since (sI+A)^"*^ G H^^^ and H(s) G 7i^^^ , it follows from Lemma [2] that the spectral radius 
p{s) in (fTSll is a subharmonic function on C+. As a direct consequence, we have, among all, that 
p{s) is a continuous bounded function on C+ and satisfies the maximum modulus principle (see, e.g., 
|15j and references therein): p{s) achieves its global maximum only on the boundary of C+. Since 
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p{s) is strictly proper in C+, i.e., /?(s) — > as \s\ +00 while keeping s £ C+, it follows that 
suPsgc+ P i^) < sup^gjj p {joj) . According to the latter inequality, condition is satisfied if 

P U^) = P [{jtol+^r' H(ic^)) < 1, \/ioe M\{0}. (19) 

Denoting by || A||g^ = max^ | [AJ^,^ | the maximum row sum matrix norm and using p (A) < || Ajj^^ 
VA G C^>^^^ [18], we have 



pijuj) < (ja;I+A)-iH(ja;) 



:max < max — — — — — < 1, 



r 



where in the last inequality we used ([3D [see assumption a2)]. Since in the second inequality, the 
equality is reached if and only if w = 0, we have p (jco) < 1 for all 7^ 0, which guarantees that 
condition (jlSp is satisfied. Hence, according to Lemma [H given any set of initial conditions {(j)^} 
satisfying a4), the trajectories ^{t) as t ^ +00, with exponential rate arbitrarily close to 

r = — {miuj [ Re{sj}| : p{si) = and Sj / 0}, where p{s) is defined in (fT2|) and S TZ{1n) (because 
of L^'"" = and rank(L) = — 1). In other words, system (|14p exponentially reaches the consensus 
on the state. 

Necessity: The necessity of quasi-strong connectivity of the digraph for the network to reach a global 
consensus can be proved as in [11^ Appendix C.2] by showing that, if the digraph associated to ^ 
is not QSC, different clusters of nodes synchronize on different values [11, Corollary 1]. This local 
synchronization is in contrast with the definition of (global) synchronization, as given in Definition 1. 
Hence, if the overall network has to synchronize, the digraph associated to the system must be QSC. 

References 

[1] R. Olfati-Saber, J. A. Fax, R. M. Murray, "Consensus and Cooperation in Networked Multi-agent Systems," in Proc. of the 
IEEE, vol. 95, no. 1, pp. 215-233, Jan. 2007. 

[2] W. Ren, R. W. Beard, and E. M. Atkins, "Information Consensus in Multivehicle Cooperative Control: Collective Group 
Behavior Through Local Interaction," IEEE Control Systems Mag., vol. 27, no. 2, pp. 71-82, April 2007. 

[3] R. Olfati-Saber and R.M. Murray, "Consensus Problems in Networks of Agents with Switching Topology and Time-Delays," 
IEEE Trans, on Automatic Control, vol. 49, pp. 1520-1533, Sep., 2004. 

[4] M.G. Earl, and S.H. Strogatz, "Synchronization in Oscillator Networks with Delayed Coupling: A Stability Criterion," P/ij/sicaZ 
Rev. E, Vol. 67, pp. 1-4, 2003. 

[5] A. Papachristodoulou and A. Jadbabaie, "Synchronization in Oscillator Networks: Switching Topologies and Presence of 
Nonhomogeneous delays,"in Proc of the IEEE ECC-CDC '05, Dec. 2005. 

[6] D.J. Lee and M.W. Spong, "Agreement With Non-uniform Information Delays," in Proc. of the ACC '06, June 2006. 

[7] J. N. Tsitsiklis, D. P. Bertsekas, M. Athans, "Distributed Asynchronous Deterministic and Stochastic Gradient Optimization 
Algorithms," IEEE Trans, on Automatic Control, pp. 803-812, Sep. 1986. 

[8] D. P Bertsekas and J.N. Tsitsiklis, Parallel and Distributed Computation: Numerical Methods, Athena Scientific, 1989. 

[9] V. D. Blondel, J. M. Hendrickx, A. Olshevsky, and J. N. Tsitsiklis, "Convergence in Multiagent Coordination, Consensus, 
and Flocking," in Proc. of the IEEE CDC-ECC'05, Dec. 2005. 

[10] S. Boyd, A. Ghosh, B. Prabhakar, D. Shah, "Randomized gossip algorithms," in IEEE Trans, on Information Theory, Vol. 
52, pp. 2508-2530, June 2006. 

[11] G. Scutari, S. Barbarossa, and L. Pescosolido, "Distributed Decision Through Self-Synchronizing Sensor Networks in the 
Presence of Propagation Delays and Asymmetric Channels," in IEEE Trans, on Sign. Proc, vol. 56, no. 4, pp. 1667-1684, 
April 2008. 



13 



[12] S. Baxbarossa, and G. Scutari, "Decentralized Maximum Likelihood Estimation for Sensor Networks Composed of Self- 
synchronizing Locally Coupled Oscillators,", IEEE Trans, on Sign. Proc, vol. 55, no. 7, pp. 3456-3470, July 2007. 

[13] S. Barbarossa and G. Scutari, "Bio-inspired Sensor Network Design: Distributed Decision Through Self-synchronization," 
IEEE Sign. Proc. Magazine, vol. 24, no. 3, pp. 26-35, May 2007. 

[14] W. Rudin, Real and Complex Analysis, McGraw-Hill, International Student Ed., 1970. 

[15] S. Boyd and S. A. Desoer, "Subharmonic Functions and Performance Bounds on Linear Time-Invariant Feedback Systems," in 
IMA Jour, of Math. Control & Inf., Vol. 2, pp. 153-170, 1985. 

[16] R. Bellman and K.L. Cooke, Differential- Difference Equations, New York Academic Press, 1963. 

[17] K. Gu, V.L. Kharitonov, J. Chen, Stability of Tyme-Delay Systems, Control Engineering Series, Birkhauser, 2002. 

[18] R. A. Horn and C. R. Johnson, Matrix Analysis, Cambridge Univ. Press, 1985. 

[19] G. Mergcn, L. Tong, "Type-based Estimation over Multiaccess Channels," in IEEE Trans, on Signal Proc, Vol. 54, pp. 
613-626, Fcbr. 2006. 

[20] L. Pescosolido and S. Barbarossa "Distributed Decision in Sensor Networks based on Local Coupling through Pulse Position 
Modulated Signals," in Proc. of IEEE ICASSP 2008, March 30 - April 4, 2008, Las Vegas, NV, (USA). 

[21] L. Pescosolido, S. Barbarossa, and G. Scutari, "Radar Sensor Networks with Distributed Detection Capabilities," in Proc. of 
IEEE Radar Conference 2008, Sheraton Golf Parco dci Medici, May 26-30, 2008, Rome, Italy. 

[22] S. Barbarossa, G. Scutari, A. Swami, Achieving Consensus in Self-Organizing Wireless Sensor Networks: The Impact of 
Network Topology on Energy Consumption, in Proc. of IEEE ICASSP 2007, April 15-20, Honolulu, Hawaii (USA). 



14 



